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Abstract 

The effect of the random quantum transverse field VL on the tricritical behavior of the 
spin-1 Blume- Emery- Griffiths (BEG) model is studied using effective field theory. 
It is found, that the tricritical behavior depends on both the biquadratic interaction 
K., single- ion anisotropy A, and the concentration p of the disorder of Vl. Indeed, 
there exist a special value pi of the probability p below which the tricritical behavior 
disappears. In addition, at sufficiently low temperature, the system exhibits long-range 
order accompanied by the tricritical behavior below a special value p2 of the probability 
p. 
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1- Introduction 

The Blume- Emery- Griffiths (BEG) model [1] is a generaUzation of the spin- 1 Ising 
model. It considers both bilinear (J) and biquadratic (K) nearest- neighbor pair in- 
teractions, in which a single- ion anisotropy parameter A is included. With vanishing 
bi-quadratic interactions, the model is known as the Blume-Capel model [2,3]. The 
BEG model has attracted a great deal of attention since it was originally proposed 
to describe phase separation and superfluid ordering in ^He — ^ He mixtures. Since 
its introduction, the BEG model has been extended using a variety of techniques, to 
solid-liquid-gas systems, multi-component fluid liquid crystal mixtures [4,5], magnetic 
materials [6 — 8], critical behavior and multi-critical phase diagrams [9 — 12], the re- 
entrant phenomenon [13 — 20], and to study metastable and unstable states [21 — 23]. 
Furthermore it is the simplest model that can be used for modelling the behavior of the 
liquid, solid and vapour phases of a real materials [24, 25]. The phase diagrams of the 
BEG model for < K/ J have been studied by the mean field approximation [26 — 28], 
the position space renormalization group method [29], the cluster variation method 
[30 — 31], series- expansion methods [32], the transfer- matrix method [33], the con- 
stant coupling approximation [34], linear- chain approximation [35], and on the Bethe 
lattice using exact recurring equations [36]. On the other hand the BEG model with 
repulsive biquadratic coupling i.e., K/ J < was a subject of interest of many authors 
[37 — 41], these studies have shown a variety of interesting features as for example single 
and double reentrancy regions and ferrimagnetic phases. Keskin and Ekiz [42] have 
also investigated the thermal variations of the BEG model with repulsive biquadratic 
interactions by using the lowest approximation of the cluster- variation method. 
The quantum effects on the tri critical behavior of the spin-1 BEG model are very im- 
portant. However, many works have been done in this context [43 — 48]. Particularly, 
the effect of random fields on the phase transitions in quantum Ising systems has been 
a subject of much interest [48 — 56] , since they show various types of multicritical phe- 
nomena. 

Our aim in this paper, is to study the effect of the random quantum transverse field 
on the behavior of phase diagrams of the spin-1 BEG model using the finite cluster 
approximation [57 — 59] . It is found, that the tricritical behavior disappears below a 
concentration pi of the disorder of Q. In addition, at sufficiently low temperature, it is 
shown that the system exhibits for any value of Q, a long-range order accompanied by 
the tricritical behavior, below a special value of the probability p2- Pi and p2 depend 
on both biquadratic interaction K and single site anisotropy A. 

The paper is organized as follows: in section 2 the model and the method are well 
described; in section 3 results and discussions are given; section 4 is reserved to the 
conclusion. 
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2- Model and method 

The spin-1 BEG model in a transverse field is described by the following Hamiltonian 



(1) 



Where J and K represent respectively the bilinear and bi-quadratic couplings, A is 
the single site anisotropy, indicates summation over nearest-neighbor sites of the 

simple cubic lattice, J2i indicates summation over all the sites of the lattice, while Si^ 
and Siz are the x and z-components of the spin-1 Pauli matrix respectively given by: 




fli is the transverse magnetic field, which is given by the following distribution law 

p(a) = p6{n, + p)6{n,) (2) 

Where p is the transverse field concentration and S{x) is the Dirac function defined as 
follows: 

5(0) = l;5(x^0) = (3) 

Using a single cluster approximation method in which attention is focused on a 
cluster comprising just a single selected spin labelled 0. And the neighboring spins 
with which it directly interacts, then the Hamiltonian containing namely Ho is given 

by 



Ho^ASoz + {A + C)Sl-noS, 



Ox 



Where 



- J0;C ^ -Kt] 



(4) 
(5) 
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With 

(6) 

i=i j=i 

are operators functions of Sj^ {j = 1, ...,N); N is the coordination number. 
This single site Hamiltonian can readily be diagonalized and its eigenvalues and eigen- 
vectors are easily obtained. The three eigenvectors corresponding to the eigenvalues 



With 



And 



Are 
With 



Xk ^'^{A + C+^cos{ifik)y,k^ 1,2,3 (7) 



(f^ = ^arccos{-27q/2g) + 2/3{k - l)7r (8) 



r = + _ (A + C)V3 (10) 

g = -(A + C)(2A'-^(A + C)'-Q')/3 (11) 

9 



>= Q;fc|+ > > +7fe|0 > (12) 



^ |1^(A,-A-C + A)| 

2Vyi]2((^^ _ A - C)2 + A2) + ((Afe - A - C)2 - ^2)2 

7fc = V2 Qffc (15) 

In a representation in which S^^^ is diagonal. The starting point of single-site cluster 
approximation is a set of formal identities of the type 

rro(go,)'exp(-/;/fo) 
■= '^O") Tr„e.p(-/?/f„) 
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Where {SoaY is the a- component of the spin operator Sq raised to the power 1, 
< ('S'oa)' >c denotes the mean value of {SoaY for a given configuration c of all other 
spins, i.e when all other spins, Si{i 7^ 0) , have fixed values. Ttq means the trace 
performed over 6*00 only. /? = 1/KbT; T is the absolute temperature and Kb is the 
Boltzmann constant. 



The equations (16) are not exact for an Ising system in a transverse field, they 
have nevertheless, been accepted as a reasonable starting point in many studies of that 
system [60]. 

To calculate < Soa >c and < (^oa)^ >c, one has to effect the inner traces in eqs. (16) 
over the states of the spin and this is the most easily performed using the eigenstates 
of equation (12). By setting I — 1 and I — 2 in eq.(16) it turns that 



E LiK-/3k)exp(-/3 AO 
ELiexp(-/3Au) 



<>^Oz>c— — ^33;^ — ^v^; U'j 



^ q ^ . AT Ek=i(«k + /3k)exp(-/3Ak) . . 

< ^Ox >c = Vi ^^3 (18j 

Ek=i exp(-/3Ak) 

_ Ek=i(«k + /gk)exp(-/jAk) . . 

Ek=i exp(-/3Ak) 



ELi((«k_+A£(2 + 7|)exp(-/3Ak) 
ELiexp(-/3Ak) 

The magnetizations ma{cK — z,x) and the quadrupolar moments ga(Q; = z,x) are 
given respectively by 



spins 

(21) 

Qa =<< 9a{0: V, ^o) > D> spins (22) 

With 

fa{0,V,^o)^<Soa>c (23) 



ga{0,vM^<{Soaf>c (24) 

Where < ... > spins means the average over all configurations of the nearest neighbors 
spins Si{i 7^ 0) of the site , and < ... >d denotes the average over all configurations 
of the disorder of the transverse field. Using the distribution of Qi mentioned above, 
the average over the disorder of fli of and are given by 
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(25) 



< ga{9, T], no) >D^ j 9a{e, Tj, Oo)P(^^o)ci^^o (26) 

However to calculate and we use an expansion technique of the spin- 1 cluster 
identities [59 ]. 

Suppose one considers the general product Ilili J2t=o ^^^^ contains 3^ terms. From 
these terms one may collect together all those terms containing u factors of Sf^ and v 
factors of 5'^^. Such a group is denoted by {S"^, Sz}n,u,v For example, if N = 3, u — 1 
and V — 2, then 



{Si, Sz}3S^2 — Sf^S2zSsz + S2zSizS3z + Sl^SizS2z (27) 

Our aim now, is to expand /a(6', ^, f^o) and gaiO, ^, ^^o) in terms of these {Si, Sz}n,u,v 
Thus if one writes 

N N-v 

U9, T], llo) = E E AS,v(N, ^lo){S^ S}n,u,v (28) 

v=0 u=0 

N N-v 

ga{9, ri, Oo) = E E B^,v(N, ^^o){S^ S}n,u,v (29) 

v=0 u=0 

The problem is to find the coefficients A^.^{N, Omegao) and B'^.^{N,Omegao). To 
achieve this is advantageous to transform the spin-1 system to a spin- 1/2 representation 
containing the Pauli operators Ciz — 1, —1. This may be accomplished by setting 
Siz — Tiz(7iz with Tiz = 0, 1. In this representation, fa and ga become 

N N-v 

Ue, V, l^o) = E E Au,v(N, ^^o){tz, T,ajN,u,v (30) 

v=0 u=0 
N N-v 

ga{e,v,^o) = E E BS,,(N,Qo){t„T,(7Jn,u,v (31) 

v=0 u=0 

And must hold for arbitrary choices of Tiz- Suppose one now chooses the first r out of 
the N operators r^z to be unity, and the remainder zero. Then (30) and (31) become 

fa{e,v,no) = EEA°v(N,i^o)cr{^4r,v (32) 

v=0 u=0 
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go,{e, ri, Qo) = E E BS,v(N, Jlo)Cr {^z}r,v (33) 

v=0 u=0 

Where {(Tz}r,v is the sum of all possible products of v spin operators, aiz, out of a 
maximum of r, and the C™ are the binomial coefficients m!/n!(m — n)\. That's, 

U0, T), Qo) = E K (r, ^0)W.}r,. (34) 
v=0 

ga(e,ri,no) = E<(r,^^o){^z}r,v (35) 

v=0 

With 

6:(r,Oo) = EAS,,(N,Oo)Cr (36) 

u=0 

<(r, Qo) = E B^,v(N, no)Cr (37) 

u=0 

The spin-1 problem of eqs. (28) and (29) containing N spins has thus been transformed 
to a spin-1/2 problem containing r spins [59]. The advantage of doing this is that now 
enables one to use directly the results already established in [61] for the spin 1/2 
system. Specializing the results of [61] to a single group of r spins, one has for the 
current problem 

6^(r, Oo) = l/(2^Cr/) E C[ei(r, v)fi«(r, Qq) (38) 

i=0 

<(r,0o) = l/(2'-C/)EC[ei(r,v)giJr,no) (39) 



i=0 



Where 



And 



e^{r,v) = Y.{-^r C^Cl-J, (40) 

M=0 



/,„(r,(]o) = /a(r-2?:,Qo) (41) 
gia{r,flo) ^ ga{r -2i,no) (42) 

Once the coefficients 6"(r, (^o) ^-nd (i"(r, (^o) have been calculated, the coefficients 
^m,d(A^) ^o) and B",^j{N, VIq) may be found by the following procedure. First, Aq^^{N, VLq) 
and Bq^jj{N, Qq) are got by setting r = t> in (36) and (37). That's, 

A^^{N,no)-^K{v,no) (43) 
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And 

Bl{N,n,)^d^{v,no) (44) 

Then, the other coefficients A^^^{N,Qq) and B^,^{N,Qo) may be obtained by ex- 
pressing (38) and (39) as a reccurence relation, namely 

A'^-vA^o) = Kir, no) - AS,,(Oo)Cr (45) 

u=0 

B?-.,M^S:^\r,no)-~j:'Bl^ino)Cr (46) 

u=0 

Then the magnetizations and the quadrupolar moments Qa, where a = z,x, for an 
arbitray coordination number A^, are given by 

N N-v 

= E E < KA^o) >D< {s^, s Jn,u,v > (47) 

v=0 u=0 

N N-v 

= E E < Buv(^^o) >D< {S^ S Jn,u,v > (48) 

v=0 u=0 

Using the simplest approximation of the Zernike decoupling of the type 

< Si^Sj^ Skz---- >=< >< Sj^ > < Skz > ■■■ ioYi^ j^k^ ... 

And seeing that the number of elements of the group {S^ , Sz}n,u,v is equal to Cj^ C^^~'^ , 
then eqs. (47) and (48) become 



N N-v 



(^) = E E < Auv(^o) >D m/q,"C/C/-" (49) 



v=0 u=0 



N N-v 

<lam = E E < KA^o) >D m/q,"Cu^Cv^-" (50) 

v=0 u=0 



With 



< Al(f]o) >D= pA'^^in) + (1 -p)Al(O) (51) 

< i?:,(fio) >D= pBZ,m + (1 - p)BZM (52) 

A first order transition is characterized by the gap of the longitudinal magnetization 
ruz at the transition temperature. If the magnetization vanishes continuously at the 
transition temperature, this is the second order transition. 

3- Results and discussions 
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In the following, the ordered phase is characterized by a non zero longitudinal mag- 
netization rriz, while the disordered phase is characterized by = 0. We denote here 
after by T, Q, A and K the reduced parameters T/ J, Q/J, A/ J and K/ J respectively. 
At finite temperature, eqs. (21), (22) are solved numerically for N = 6 (the simple 
cubic lattice) using iterative method. By varying the value of the biquadratic coupling 
K, we have distinguished three topologies of phase diagrams in the (Q — T) plane for a 
fixed value of A = 3.5, which are, the existence of tricritical behavior, only the second 
order transition exists, and only first order transition is present. Hence, in the case 
K = 0.5, Fig. la shows the tricritical behavior which depends also on the probability 
of the disorder of the transverse field, that's, it exists a critical value pi below which 
the tricritical behavior disappears, this value depends obviously on the biquadratic 
coupling. Fig. lb shows the phase diagram (Q — T), where the tricritical behavior per- 
sists. So We can say that by fixing the biquadratic coupling value K, the probability 
may determine the transition order, as shown in the line corresponding to p = 0.2 in 
Fig. la, in which just the second order transition appears, while there is no first order 
transition. 

One feature should attract our attention from these two figures, there is the fact that 
the first order transition is more dominant in the Fig. lb. In principle, we can say that 
decreasing biquadratic coupling favors the first order transition. In fact. Fig. Ic shows 
the phase diagram (f2 — T) for K — 0.2, in which only the first order transition appears 
for any probabihty p. In contrast, for X = 3 as is seen in Fig. Id, the (Q — T) gives 
rise to just the second order transition. More precisely, there exist a special value of 
the biquadratic interaction Ki (= 0.27 for p = 0.4) below which there is no tricrit- 
ical behavior (only first order transition can take place), and when K increases Tfri 
decreases, up to a value oi K — K2 {— 0.72 for p — 0.4) above which the tricritical 
behavior disappears (only second order transition can give occurrence), see Fig. 2a. 
In contrast, by fixing K and increasing p, Tt^i increases, see Fig. 2b. Such that, there 
exist a special value pi ( = 2.252 for K = 0.5) of the concentration of the disorder of 
fl, below which the tricritical behavior disappears. 

The behavior of the tricritical transverse field as a function of K and p is given in Figs. 
3a and 3b respectively. It is found that D,tri increases when increasing K (Fig. 3a), 
while it decreases with p as shown in Fig. 3b. 

Besides, for any value of Q and at sufficiently low temperature, the system exhibits a 
long-range order accompanied with a tricritical behavior below a special value of the 
probability p2- 

4- conclusion 
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In this work we have investigated the effect of the random transverse field f2 on the 
tricritical behavior of the quantum spin-1 Blume- Emery- Griffiths (BEG) model. We 
have found in the {Q, T) phase diagrams, that, the tricritical behavior depends on both 
the biquadratic interaction K, and the concentration of the disorder p of fl. It is also 
found, that, at sufficiently low temperatures, the phase diagrams exhibit long-range 
order accompanied by the tricritical behavior. 
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Figure captions 



Fig.l Phase diagrams in the {^,T) plane for A = 3.5 

(a) : 

(b) : 

(c) : K = 0.2 

(d) : K = 3. 

Fig.2 (a): The behavior of the tricritical temperature as a function of K for p = 
(b) : The behavior of the tricritical temperature as a function of p for K — 

Fig.3 (a): The behavior of the tricritical transverse field as a function of K for p 
(b) : The behavior of the tricritical transverse field as a function of p for K 
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Fig. 1c 



■ First order transition 
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Fig. 1a 
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